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Abstract 

By using variational calculus and exterior derivative formalism, we 
proposed in [30] and [11] a new geometric approach to electromagnetism 
in pseudo-Finsler spaces. In the present paper, we provide more details, 
especially regarding generalized currents, the domain of integration and 
gauge invariance. Also, for flat pseudo-Finsler spaces, we define a gen- 
eralized energy-momentum tensor consisting of two blocks, as the sym- 
metrized Noether current corresponding to the invariance of the field La- 
grangian with respect to spacetime translations. In curved spaces, one 
of the blocks of the generalized energy-momentum tensor is obtained by 
varying the field Lagrangian with respect to the metric tensor and the 
other one, by varying the same Lagrangian with respect to the nonlinear 
connection. 



1 Introduction 

Classical electromagnetism is one of the most "rounded" theories of physics 
and one has strong reasons to say that there is little to be added to it. There 
exist several beautiful geometrical descriptions of this theory in classical general 
relativity. 

Still, what if spacetime is described not by Riemannian, but by Finslerian 
geometry? As G. Yu Bogoslovski and H. Goenner noticed, "spacetime may be 
not only in a state which is described by Riemann geometry but also in states 
which are described by Finsler geometry" . A lot of authors have already consid- 
ered Finslerian models for spacetime (see, for instance, [1], [l]-[2], [8]- [10], [15] . 
[H], [10], [12], [13], [H], [H], [33], in order to cite just a few of them). 

Regarding electromagnetism, we can expect that in spaces with Finslerian 
geometry, the corresponding equations would change and that we might even 
have to deal with some new quantities. In this paper, we will investigate from 
a mathematical point of view these possible changes. 
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The first idea we must have in mind is that in spaces with Finsler geometry, 
the metric tensor depends on the directional variables. Since Maxwell equations 
involve the metric tensor, we notice that the solutions (hence, the electromag- 
netic tensor) may also depend on these. This means that the natural space to 
work on is not the spacetime manifold M, but its tangent bundle TM. 

Thus, we extend the classical ideas of electromagnetic field theory to Finsler 
spaces, as an application of the geometry of the tangent bundle TM . 

Though we will speak throughout the chapter about Finsler spaces, all the 
theory remains valid, with minimal changes, for more general anisotropic spaces 
(Lagrange, generalized Lagrange ones). 

This paper is a continuation of [TT] and [3D]. It is based on variational 
calculus and classical methods in theoretical physics (adapted to the tangent 
bundle). This approach offers an alternative to the existing one by R. Miron 
and coUaborators, [T7], [H], 

2 A brief overview of the Riemannian case 

There are multiple definitions of the electromagnetic tensor, points of view and 
formulations of the basic equations of electromagnetism on Riemannian man- 
ifolds. Namely, the electromagnetic tensor can be regarded as the curvature 
of a line bundle over the given manifold or it can be described in terms of 
nonlinear /linear on TM as in [TS], or in terms of differential forms. 

In the following, we will adopt the language of differential forms, as it is the 
most tightly related to variational calculus. 

Let us consider a pseudo- Riemannian manifold {M, g) of dimension 4, thought 
of as spacetime manifold. We denote local coordinates on M by x = {x^)i^Q-3 
and use the numbering from to 3. The first coordinate is regarded as the time 
coordinate and {x°')^^y3^ ^ spatial coordinates. As required by general rela- 
tivity, the metric g = g{x) is supposed to have Lorentz signature (-1-,—,—,—). 
Here are some other notations and conventions we will use in the following: 

- Latin indices i, j, k, ... take values from to 3; Greek indices a, /3, 7, ... take 
values from 1 to 3; 

- for a vector field v — (w*)j^Q-3 on M, v will denote the spatial vector 

V = (^^")„=o- 

. , , . . . , d 

- k ~ partial derivative with respect to t^— r; 

ox'^ 

d - 

- \k ~ Levi-Civita covariant derivative with respect to -t—t; 7* ,i. - Christoffel 

ox'^ ■' 

symbols of 

- g = det ((jij); whenever it is not clear from the context whether we refer to 
g as the metric tensor or to the determinant of the corresponding matrix, we 
will specify this; 

- b : TM — ^ T*Af, jj : T*M — >• TM - musical isomorphisms (lowering/raising 
indices); 

- d - exterior derivative of differential forms, * - Hodge dual; 
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- - signature of the permutation (ii, ip); 

- V - gradient taken with respect to the spatial coordinates (x"); 

- d'^x = dx° A dx^ A dx"^ A dx^, d^x = dx^ A dx^ A dx^; 

- dil — \/\g\d'^x, dV — ^^^^ d'^a; - the invariant Riemannian volume element 

V5oo 

on spacetime and on the spatial manifold respectively. 



2.1 Distances, volumes, divergence, codifFerential 

Let us remind for the beginning some very quick facts about computation of time 
intervals, distances and spatial volumes in general relativity (|16], pp. 315-320). 

The (squared) arclength element ds^ on the spacetime manifold M can be 
written as 

ds'^ = gooidx")'^ + 2goadx°dx" + gapdx°'dx^ . 
The spatial arclength element is defined as 

dl =^^pdx dx'\ ^„^-g^p-\ , a,/^ = l,3. 

500 



The determinant of the spacetime metric g is 

9 = -5007; 

we only consider reference frames for which both the determinant of the spatial 
metric 7 and goo are positive: 
7 := det(7,y ) > and goo > 0. 

Here are some other relations we will use in the following. 
The partial derivatives of y^jgj are given by: 



diln^\g\) 1 



dx^ 2 



= 0.9 9ih,j = 7 



Id -I 

Consequently, the divergence diviV") — — =— ^(F*v|g|) of a vector field can 
be written in terms of covariant derivatives, as 

div{V) = V\. 

Also, by expressmg as — f-^— 7, we get 

d{\^^\) = -\g,,dg'^K (1) 

The latter equality is particularly useful when varying Lagrangians with 
respect to the metric. 
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The codifferential of 5£, = (—1)^ * of a p-form ^ = ~f^iii2...ip^^'^ ^ ■•• ^ 

dx^p is a (p — l)-form, with the property {ri,d^) — {drj,^) , where ( , ) denotes 
the inner product of p-fornif0. For a 2-forni, it is given by 

2.2 4-potential and electromagnetic tensor 

The 4-potential is geometrically described in classical general relativity as a 
1-form 

A^A,{x)dx\ (2) 
The electromagnetic tensor (or Faraday 2-form) is described as the 2-form 

F = dA. (3) 

In local coordinates, this is 

F ii^,7c^ia:^ A c^a:^ (4) 

where 

Fjk = Aky ~ Aj\k- (5) 

Due to the symmetry of the Levi-Civita connection, the latter expression 
can be actually written in terms of partial derivatives only: 

F,k = ^fc,, - A,,k. (6) 

In the language of differential forms, the homogeneous Maxwell equations 

Fij\k + Fki\j + Fjk\i — 0. (7) 

become 

dF = 0. (8) 

Remark. There exist two possible approaches regarding the potential A 
and the electromagnetic tensor F. 

1. One can consider as a fundamental object the electromagnetic tensor F^ 
regarded as a closed 2-form on the manifold. In this case, the homogeneous 
Maxwell equation dF — , i.e., the closure condition for F, is taken as 
an axiom. If the manifold M is topologically "nice enough", then one 
can apply Poincare's lemma, which entails the existence of a 1-form A, 
such that F = dA. That is, the existence of the 4-potential is seen as a 
consequence of the homogeneous Maxwell equations. 



^The inner product of two p-forms d = Oi-^...ip^^^ A ... A e'p and i/i = '/'ji...jp6"'^ A ... A e^p 
is given by J g*!-*! ■■■g^P-'P 9ii...ip'4>ji jp'^^' where the integral is taken on the whole manifold 
(assuming that the integrands have compact support). 
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2. Some authors consider the potential l-form A as a fundamental object and 
define F as its exterior differential; in this approach, the homogeneous 
Maxwell equation is obtained as an identity. 

Actually, for a Lagrangian theory of electromagnetism, it is essential to have 
both a l-form A and a 2-form F, related by ([8]). 

An important property of the electromagnetic field is gauge invariance. 
Namely, the field strength tensor F is invariant to transformations 

where V' : -A/ — ^ M is a differentiable function. 



2.3 Lagrangian, equations of motion and inhomogeneous 
Maxwell equations 

The second pair of Maxwell equations (inhomogeneous Maxwell equations) and 
also, the equations of motion of charged particles in a given electromagnetic 
field are obtained by variational methods. 

The total action attached to the field and to a system of particles is 

S = -^mc j ds -^'^ j Ak{x)dx^ - j F.jF'^dn. (9) 

Here, m, q, c are constants (m denotes the mass of a particle, q, its charge, c, 
the speed of light in vacuum), dVl — yj\g\d'^x is the invariant volume element on 
spacetime and the sums are taken over the particles in the system. The volume 
integral is taken over a bounded interval of time and over the whole spatial 
manifold, under the assumption that far away from sources, the field vanishes. 
Thus, we can actually think the integral as taken over a " large enough" compact 
domain in M . 

The first term, Sp := — ^mc / ds, corresponds to the Lagrangian Lp := 

— ^ mods of the free particles. 

The second term Sint '■— ^'12~ I ^k{x)dx'^, given by the Lagrangian Lint ■= 

— ^ -Ak{x)dx'' characterizes the interaction between the particles and the field. 

The third term Sf :— — — f FaF'^^dil — — f F* F d'^x, characterizes the 

■' , IGttc . Wire 
electromagnetic field in the given curved space. 

By keeping the electromagnetic field fixed and varying the trajectory of a par- 
ticle in the action S (which actually means varying the trajectory in Sp -\- Sint), 
one obtains the equations of motion of particles subject to both gravitational and 
electromagnetic field, i.e., the expression of the Lorentz force. If, conversely, in 
the action S we keep trajectories fixed and vary the electromagnetic field (which 
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actually means to vary the electromagnetic field in the sum Smt + Sf), we get 
the field equations, i.e., the second pair of Maxwell equations. 



Let us notice that Sp and Sint are line integrals, while Sj is given by a 
volume integral. Hence, if we want to vary Smt + S"/, we have to write this 
sum as a single volume integral, too. This is achieved by means of the notion 
of charge density. 

Charge density p is defined as the amount of electric charge in a given spatial 
volume and it is basically a function of time and spatial coordinates: 

P = P{x)- 

The integral of p over a certain region of space provides the total charge 
situated inside that region: 

q = I pdV, (10) 



where dV = — d x is the spatial volume element. In this writing under an 

integral, it is supposed that we actually see the charge distribution as "contin- 
uous" . Total charge is invariant to coordinate changes. 

For a discrete distribution of charges gi, ...,(/„ in a given volume, we can still 
write the total charge in the form of the integral (|10p . if we define the charge 
density by means of the Dirac delta function: 

n 

where x = {x'^,x'^,x^) and X(„-) is the position vector of the charge qa- 
By using relation (|10p . Sint is written as 

S^nt = I AJ'dn, 



where the quantities 



J' = (11) 



are the components of a vector field J, called the 4-current. 

Thus, the sum Si := Sint + Sf can be written as a single integral as 



Si = - j il^r + -l-F,^F^i)^\d^x. 



By varying the above Lagrangian with respect to the potential A, one gets 
the field equations, i.e., the inhomogeneous Maxwell equation: 

SF = - — J,, (12) 

c 
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or, in local writing, 



Remark. The Maxwell equations and eventual choices of ip in the transfor- 
mations A A -\- dip do not completely determine the potential A. Hence, to 
A, one can still impose supplementary conditions (gauges). The most common 
is the Lorenz gauge ^*|. = 0; under this condition the inhomogeneous Maxwell 
equations become: 

, Ajm _ 4^ 



A'^'.+A^R% = - — J\ (14) 



where R^j = g^^R^jk components of the Ricci tensor. 

By using inhomogeneous Maxwell equation, one obtains that the 4-current 
J identically satisfies the continuity equation: 

div{J) = d{*J\,) = 0. (15) 

i.e., 

J\i = 0- (16) 
From a physical point of view, the continuity equation is equivalent to the 

charge conservation law. 

Let us now consider a single particle, subject to the action of a given (fixed) 
electromagnetic field and determine the trajectory of this particle. This can be 
achieved by varying the action S with respect to the trajectory. That is, we 
actually have to vary the action ^2 := Sp + Sint, written as a line intcrgral. 

We notice that that the integral Sp + Sint does not depend on the choice of 
the parameter on the path of integration. Thus, we can choose this parameter 
according to our wish. So, let us choose the arclength s as a parameter. With 
dx 
ds 



this choice, we have 
Thus, we can write 



/ ds 

= VdijX^X^ = — = 1. 



^2 = — mc J ds — - J Ai{x)dx^ = — J {mC'^gijX''X^ + AiX^)ds, (17) 

where the integral is taken on some fixed compact interval [so, si]. 
The Euler-Lagrange equations for the above Lagrangian are 

= 1f^xJ, i = o,3, (18) 



Dx'' dx^ 

where — — = — — h j^^.x^x'' is the Levi-Civita covariant derivative. 



ds ds 
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The right hand sides of the above equations provide the expression of the 
Lorentz force in the given curved space. Also, the first two terms in ^ provide 
the canonical momentum 

Pi = wct^ + -Ai. (19) 

||x|| c 

2.4 Energy-momentum tensor 

Another important quantity in general relativity, is energy-momentum tensor 
(or stress- energy tensor) T. In classical general relativity, the energy-momentum 
tensor is symmetric and, for a closed system, its covariant divergence is zero. 

A. In fiat Minkowski space 

In the Minkowski space (M^, rf) {rj = diag{l, —1, —1, —1)) of special relativity, 
it makes sense to speak about spacetime translations 

X I— > a; + a (a - constant 4- vector). 

The Lagrangians L/, Lint, Lp in (|9]) are all invariant with respect to these trans- 
lations (this can be easily checked, by noticing that neither of them depends 
explicitely on the spacetime coordinates x^). 

Moreover, in this case we have g = det{ri^j) = — 1 and covariant derivatives 
coincide with partial ones. 

According to Noether's theorem, the invariance of an action 



to translations implies that the quantities 

are conserved (divT — 0). They define a tensor of rank two (the Noether current 
attached to the Lagrangian). 

The Noether current is generally not symmetric. Still, this situation can be 
"mended" by adding a divergence term: 



(where ip^''\x) — —ip^''''{x) are functions of class at least two), which does not 
affect the value of the action integral (assuming, as usually, that on the boundary 
of the integration domain, the involved functions vanish). Thus, one obtains a 
symmetric tensor T of rank two, with 

dx'^ 
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The energy-momentum tensor of the electromagnetic field in flat space is 
defined as the symmetrized Noether current given by the invariance of the field 
Lagrangian Lf to spacetime translations. 

For electromagnetism, we have = j4(i.) and 

167r ^ 

By supposing, at first, that p = (which implies J = 0), we have: T\ = 

_L(_i^'feAfe., + ^SlF^kF^'"); by adding the quantity -^A' = -^-(^^f 
Att 4 47r ' 4tt ' 

one gets the energy-momentum tensor as: 

T\ = ^{-F"'Fik + \s\FjkF^''). (20) 

Thus, if J = 0, then 

div{T) = 0. 

In the situation when we have charged matter ( J ^ 0) , the energy-momentum 
tensor satisfies the identities: 

= --F^jj' (21) 

(which can be proved by means of Maxwell equations). In brief, 

div{T) = -UjF. (22) 

The quantity -ijF is called the density of Lorentz force, 
c 

B. In curved spaces 

In general Relativity, the energy-momentum tensor T is defined by the rela- 
tion 

^9S=yJ TikSg^'dQ = -^ J T^'^Sgikdn. (23) 
Variation with respect to the metric leads to: 

Tij = ^ i-Fj "Fik + \gijFikF'% (24) 
which agrees to the expression of the energy-momentum tensor in flat space. 

Again, by using (both homogeneous and inhomogeneous) Maxwell equations, 

one gets that in curved pscudo-Ricmannian spaces, the covariant divergence of 
the stress-energy tensor of the electromagnetic field is equal to minus the density 
of Lorentz force: 

^'i\J=-\P^^j'- (25) 

Conclusion. In a geometric language, the above fundamental equations of 
electromagnetic field theory can be written briefly as: 
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• F = dA. 

• dF = - homogeneous Maxwell equation; 

An 

• SF — J\, - inhomogeneous Maxwell equation; 

c 

• div{J) = - continuity equation; 

• V • r = — ijF - energy-momentum conservation (where V • T denotes 

c 

covariant divergence) . 



3 Some geometric structures in Finsler spaces 

Let, again, M be a 4-dimensional differentiable manifold of class C°°, thought of 
as spacetime manifold. This time we will also speak about the tangent bundle 
{TM, TT, M) and denote 2/*)i=(r3 the coordinates in a local chart on TM; we 
preserve the notations in the previous section, with the only difference that in- 
stead of Levi-Civita covariant derivatives, we will use other covariant derivation 
laws. Also, we denote partial derivation with respect to by a dot: .i. We will 
sometimes call the base coordinates a;* positional variables and the fiber ones, 
directional variables. 

A Finsler fundamental function on M, is a function T : TM — > R with the 
properties, [25] : 

1. F = F{x, y) is smooth for y 7^ 0; 

2. F is positive homogeneous of degree 1, i.e., F{x,Xy) = XF{x,y) for all 
A > 0; 

3. The Finslerian metric tensor: 

1 d'^F'^ 
^-^"'^^=29^' 

is nondegenerate: dei{gij{x,y)) 7^ 0, Vx e A/, y £ Tj.M\{0}. 
In the following, we will consider that the metric has signature (+, — , — , — ). 

Remark. Strictly speaking, it would be more rigorous to preserve the term 
"Finslerian" for the case when the metric tensor g is positive definite and to 
call " pseudo-Finsler" the spaces in which the metric is only nondegenerate (of 
constant signature). But since a lot of authors already use in the latter case the 
term Finsler, we will also adopt this more relaxed terminology. 



10 



In a Finsler space, the element of arc length along a curve 1 1-> x{t) is 

^, dx 
as = J^lx, — )dt. 
^ dt' 

Finsler spaces are a generalization of pseudo-Riemannian manifolds, in which 
the coefficients of the metric tensor are no longer functions defined on M, but 
on the tangent bundle TM. Actually, if on a pseudo-Riemannian manifold, the 
tangent space at each point carries a pseudo-Euclidean metric structure, in a 
Finsler space, at each fixed point x", the "norm" T{xo,y) is generally, not 
given by a quadratic form. 1-homogeneity of in y insures that the integral 
/ ds does not depend on eventual changes of the parameter along the curve 
(hence, the notion of arclength s is uniquely defined, no matter from which 
initial parametrer t we start). 

Given a (pseudo-)Finslerian metric tensor ~ gij{x,y), the corresponding 
spatial metric is defined similarly to the pseudo-Riemannian case: 7^,^ — —gai3 + 



ffOQ30/3_^ a,/3 G {1,2,3} and its determinant is det(7Q,^) — ^ ^ 



900 v9oo 
With respect to coordinate changes on the tangent bundle TM induced by 
coordinate changes (2;*) 1— (i') on the base manifold M, i.e., under coordinate 
changes 

F = x^ix), f = —^v^ (27) 

the quantities — — have a tensorial rule of transformation: — — — — — — — - , while 

oy^ oy^ ox^ oyi 

. . d 

the quantities ^— ■ have a more complicated transformation law, |12| . |15| . 

If we want to work with tensorial blocks only, then we have to use Ehresmann 
(nonlinear) connections on TM. 

We will denote by {N-'j) the coefficients of an Ehresmann connection, and 

by 

'~dx^ ^dt' '~ dt' (28) 
{dx\5y'' = df + N"'-dx^) 

the elements of the corresponding adapted basis and of its dual cobasis. Thus, 
with respect to coordinate changes (P7|) . Si and 9j have tensorial rules of trans- 
dx^ ^ g dx^ 
dx^ ' * dx^ 



formation, i.e., 5i — ^ . 5ji, di = _ dji . 



In the adapted basis, any vector field V on TA4 can be written as V 
V^Si + V^di] the component 

hV = V'5i 

is a vector field, called the horizontal component of V, while 

vV = V'd, 
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is also a vector field, called its vertical component. Similarly, a 1-form cj on 
TM can be decomposed into invariant blocks as a; = ujidx^ + WtSy^, with 

hu = ujidx^ 

called the horizontal component, and 

vuj — ujfSy^ 

the vertical one, |15) . Accordingly, any tensor field on TM is decomposed with 
respect to the Ehresmann connection into invariant blocks. 

In the following, whenever needed to make a clear distinction, we will denote 
by i,j,k, ... indices corresponding to horizontal geometric objects and by i, j, 
k, ... (with bars), indices corresponding to vertical ones - though, unless needed, 
we will not be too strict in this respect. By capital letters A, B,C, ... we will 
always denote indices which take values corresponding to both distributions: 
A,B,C,... e {i,j,k,...,t, j, k,...} 

Let us complete g up to a block metric (an hv-metric, [151) : 

GAB{x,y) ^ gij{x,y)dx' (S) dx' + Vjj{x,y)5y'- Sy' . (29) 

where o is the given Finslerian metric tensor and u is a positive definite metric 
tensoiQ. Thus, (TM, G) becomes a pseudo-Riemannian space and we can speak 
about the Riemannian (invariant) volume element on TM : 

dn ^ ^/\G\d^x A d^y. 

where G = det{GAB) (we have written d'^y instead of 6^y in the above exterior 
product, since (i^x A d^y — d'^x A 5'^y). The determinant G is, obviously, 

G ^ g - v, g = det(5ij), v = det(uy). 

The volume element dfl defines a volume element on M by: 




where = {y G T^M \ Vij{x, y)y^y^ < r^} and r = is chosen such that 

the 3-sphere of radius r in the 4-dimensional Euclidean space has the volume 
equal to 1. This volume element generalizes the idea of Holmes-Thompson 
volume in [24pl . Regarding integration with respect to x, we can assume that 

^Assuming that the topological space M is metrizable, it appears as advantageous to 
choose, for instance, a metric v which provides the topology of M. In the case when (M, g) is 
the Minkowski space, the manifold topology of M = is the Euclidean one, hence we can 
choose V as the Euclidean metric. 

''The classical idea of Holmes-Thompson volume involves integration on the indicatrix 
Ix = {y G TxM\ gijy'y-' = 1}. If the Finsler metric g is not positive definite, the indicatrix 
Ix is generally non-compact, hence this classical idea cannot be applied in our case. Choosing 
as vertical part v of the metric G on TM a positive definite one (for instance, related to the 
spacetime topology) and integrating on balls given by the metric v solves this problem. 



12 



the corresponding domain is a "large enough" compact one (in the assumption 
that far away from sources, the field is negligible and the considered time interval 
is a bounded one). 

Having a metric structure on TM, there now make sense notions such as: 
Hodge dual * or codifFerential S of p-forms on TM, gradient of a function and 
divergence of a vector field. 

The divergence of a vector field V = V^6i + V^di on TM is obtained from 
the relation CydTl = divVdQ. In the adapted frame to an arbitrary nonlinear 
conection, the divergence of a vector field is expressed as 

1 



divV = 



\G\ L 



S^{V'V\G\) + ^-^{V'V\G\) -N\.jV\ (30) 



In particular, if the vertical block of the metric is a Riemannian one v = v{x), 
then: 

- the functions 

N^-,= l\ju{x)y\ (31) 

(v) . 

where 7 j, (x) are the ChristofFcl symbols of w, are the coefficients of a nonlinear 
connection on M; 

- in terms of this nonlinear connection, the expression of the divergence is 
simplified as: 

divV = i k(l^W) + d,{V'yM)\ ■ (32) 



The codifferential of any p-form ^ = -jCiii2...ip^*^ ^ ••• ^ '^^ 

{p — l)-form S£, = (—1)^ d*; it can be also calculated from the relation 
{Vj ^0 — {d-Vj 1 where ( , ) denotes the inner product of p-formfl 
For a 2-form 

e - l^^.dx^ A dx^ + ^^Jx^ A 5y^ + ^U^V^ A 5y' 

on TM, the codifferential of ^ is a 1-forni (5^ = ijJidx^^LOa^y'^ whose contravariant 
components are given by: 



\G\ 



*The inner product of two p-forms 9 = 9i^...ij,e^^ A ... A e'p and ip = ''P A ... A e^v 
is traditionally given by f g^^-'^ ■■■g^''-'''6i-i...ip^j-^ jp*^^' where the integral is taken on the 
whole manifold (and it makes sense, for instance, for objects with compact support). In the 
case of TM, we will not integrate on the whole TM, but on a compact domain as specified 
above. 
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Choosing a nonlinear connection and a notion of covariant derivation or 
another can help to express locally in a more or less elegant form the obtained 
equations. 

It appears as convenient to choose the following linear connection DT{N), 
inspired from |15j (just - with a different covariant derivation law for vertical 
fieldsB 

= 5.X'+U,^X'\ X',^^, (33) 

- - - - r) Xi 

\t ' In ' * Qyi ' 

where 

L'jk = ^g'''iSk9hj+Sjghk-Shgjk), (34) 

The above linear connection is a distinguished connection, |15j , meaning that 
it preserves the distributions generated by the Ehresmann connection N and it 
is /i-metrical, i.e., gij\k = 0, Wjjifc = 0, '^i,j,k,i,j G {0,3}. Its only nonvanishing 
components of its torsion tensor T are 

R',k = Sy\TiSk,S,))^SkN'^~S,N\; 
P'fk - W(5,,<5,))=iVVs-^'s.- 



For the linear connection DT{N) above defined, there hold the relations: 

5, (In /M ) - L»^, , 5, (In ) " = ' (^5) 

Consequently, the divergence of a horizontal vector field — V^Si on TM 
can be written as: 

d^v{V") = {V], - P^^V^). (36) 

Another important notion for a Finsler space is the Cartan tensor C given 
by ^ 

It also has the property that 

^The choice of this linear connection instead of the classical metrical linear connection, 
|15l . appeared as a little bit more comfortable when expressing, for instance, the homogeneous 
Maxwell equation in coordinates. This is just an example. All the results can be re-expressed 
in terms of other linear connections. 
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Particular case: If w is a Riemannian metric, i.e., if 

GAB{x,y) = gij{x,y)dx^ dx^ + v^-j{x)5y^ ® Sy^ . (37) 

V 

and N is given by (PT|) . then 



and the only nonvanishing components of the torsion remain R^jf^ — ^i^jki^)y'' j 
where are the local components of the curvature of the Levi-Civita con- 
nection of V. Moreover, in this case, 

<5,(ln^) = ?V' n—O. 
This will simplify a lot of calculations. For instance, the divergence of a vector 

field V — V^Si + V^dt and the codifferential w = (5^ of a 2-form ^ = -£,i,,dx^ A 

1 2 
(ix^+^ijda;*AJ?/-' +— ^jj(5y'A(5y^ on TM are given in terms of covariant derivatives 

(IMI): 

divV = V]^ + V\ + V'C'^,, 

and 

4 4-potential 1-form 

If we want to use variational calculus in order to provide a generalization of 
electromagnetic field theory to Finsler spaces, we need a generalization of the 
notion of 4-potential. 

Let us now see how does the notion of 4-potential transform in the case when 
the geometry of the space is no longer Riemannian, but Finslerian, i.e., when 

5y = 9t]{x^y)- 

We notice that the inhomogeneous Maxwell equations involve the compo- 
nents of the metric tensor, which depend on the fiber coordinates j/Mt becomes 
clear that generally, the solutions A would depend on both x and y. Also, the 
equations themselves could become more complicated. 

Consequently, from now on, we will consider 

A^A{x,y). (38) 
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For reasons which wih be clarified later, we will also assume that the com- 
ponents Ai are 0-homogeneous in y : 

Ai{x,Xy) ^ A.i{x,y). 

That is, we will allow A to depend on the direction of y, but not on its magni- 
tude. 

In the following, we will focus on the action ^ and determine the conse- 
quences of the y-dependence of the metric g - and of the potential A. In pseudo- 
Finslerian spaces, the first term Sp formally remains the same, with the only 
difference that in the expression ds^ = gij(x,y)y^y^ dt, gij depends ony = x. 

5 Faraday 2-form and homogeneous Maxwell equa- 
tions 

Let us define the generalized Faraday 2-form (the electromagnetic tensor) in the 
same way as in Riemannian spaces: 

F = dA; (39) 

in local coordinates, this is 

F := ^Fijdx' A dx^ + Ff^dx' A V, (40) 

In terms of adapted derivatives, components of the Faraday 2-form are ex- 
pressed as 

F,j = S,Aj - S,A,, Frj=- djA, 
and in terms of covariant derivatives p3p. we get 

F,j = Aji, - F,j = -A,.-j. (41) 

In particular, if ^ = ^i^) does not depend on the directional variables, we 
get F — -{Aj^i — Aiy)dx^ A dx^ , which is similar to the expression in jTS], [TS] . 

The electromagnetic tensor F remains invariant under transformations 

A{x,y) ^ A{x,y)+dX{x), (42) 
where A : M — >■ R is a scalar function, since d{A + dX) = dA + d{dX) = dA. 

Since F is, by definition, a closed 2-form, its exterior derivative identically 
vanishes. In other words: 
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Proposition 1 There holds the generalized homogeneous Maxwell equation: 

dF = 0. (43) 

Obviously, in local coordinates, equation (|43p will acquire different forms, de- 
pending on the chosen Ehresmann connection N and linear connection DV{N). 

In terms of covariant derivatives ([M)) . equation (|43p is read as: 

Fij\k + Fki\j + Fjk\i = - ^ R'^j^F^j^; 

(«j,fc) 

Fij\k + Fki\j + Fjk-i = P%Ff.j^ - P),jFjj^, Ff,j.j + Fji,.j = 0. 

li V = v{x) and N^j{x,y) = J^jkU^'j then the second set of equations be- 
comes 

Fij\k + Fki\j + Fjk-i = 0. 

The first set in the above is the analogue (in the nonholonomic frame {Si, dj) 
on TM) of the usual homogeneous Maxwell equations. In the cases when we 
can choose the nonlinear connection N such that the horizontal distribution is 
integrable, then also the right hand sides of the first set of equations vanish. 

In the above, we have started from A as an a priori given object and defined 
F as its exterior derivative. Let us now proceed conversely and suppose that F 
is given. As we have shown in f|30|). under the assumptions that: the manifold 
M is contractible and F is a closed 2-form with vanishing Sy^ A Sy^ component, 
i.e., 

F := -Fijdx' A dx^ + Fqdx' A 5y\ dF = 0, 
there exists a horizontal form A such that F = dA. 



6 Inhomogeneous Maxwell equations 

As we have seen, in Finsler spaces the 4-potential A and the generalized Faraday 
2-form are defined on the tangent bundle TM. 

The interaction term of the total action becomes 

In the classical Riemannian case, the above integral is transformed into one 
on a domain in the spacetime M. In our case, we will transform it into an 
integral on a domain in TM. That is, we must write total charge as an integral 
(on a domain as specified above): 
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/500 

This way, Sint will be given by 



Gd-'x A d'^y. 



c J Jgoo ax^' 



With the notation 



J'-^^^ (44) 



the integral J Akdx'' is written as 

^ ^^fcda;'' = (45) 

The above expression is formally similar to the one in the pseudo-Riemannian 
case, though, here, the volume element is considered on a certain domain in the 
tangent bundle TM. 

The quantities J* (interpreted as components of the 4-current) thus define 
the horizontal component of some vector field 

J = r6^ + rOi 

on TM. 

The field equations can be obtained by varying with respect to the potential 
A the action 

Si^~- [ A — [ F * Fd^x A d^y. (46) 

c J IGttc J 

This variation leads to: 

With V = v{x) and N'^j{x, y) = 7 ^ j}^{x)y^ ^ we have 

^ {5AF''V\9\) + {F^'V\9\).j} = --J'- (48) 

c 



Notes: 1) In the integral above, in order to make sure that the expression 
has physical sense, we might need to adjust measurement units so as to have 
[Fij] — [F^j\. This can be done, by considering the fiber coordinates y* as having 
the same measurement units as the base ones (eventually, by multiplying them 
by a constant, [30]). 
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2) We remark a certain resemblance between the term {F^^^/G).-j and the 
idea of bound current in a material medium. 



Equations (j48l) gave the idea to formally generalize the inhomogeneous Maxwell 
equation as 

6F = - — J,. (49) 



c 



In local coordinates, this is: 



47r 



F'\^+r\. + Q' - - — J^ (50) 



c 

" c 



where 



Qi = -F''^P^.^+F'^aj(ln^) (51) 
= -^F^''R'^k~F'''PVk~P''Pfk- 



In particular, if w = v{x) and N^^{x,y) ~ ^7^jki^)y'^j this yields 

+ F'^^.j + F'''C\j = J* (52) 

^ \j 2 c 

The first set of equations is nothing but ([5^ obtained by variational meth- 
ods, while the second one is new. We notice the appearance of the quantities 
J' (due to both the Finslerian character of the space and the nonholonomy of 
the frame we used) which are " coupled" on TM to the usual components of the 
4-current J*. 

In the following, we will see that J' play an important role in the continuity 
equation and in the Finslerian analogue of energy-momentum conservation law. 



7 Continuity equation and gauge invariance 

Above, we have seen that 

-^J\, = 5F. (53) 

There immediately follows: SJ\, — S6F — (3, which is, diviJ) — 0. In other 

c 

words: 

^We have used the identity SSuj = {—1)"^^ {*^^ d*){*^^ d*)u} = *^^ddu} = 0. 
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Proposition 2 There holds the generalized continuity equation: 



div{J) = 0. (54) 

We have seen above that the electromagnetic tensor F is invariant under 
transformations A{x,y) i-> A{x,y) := A{x,y) + dX{x) of the 4-potential. It 
means that, in the general action the first term Sp and the third one Sf 
will also be invariant. 

The continuity equation (j54p insures that, with respect to the above trans- 
formations, Sint = — / AiJ^ y/\G\d'^x A d'^y equals Sint plus a boundary term. 

Indeed, we have (omitting the minus sign in front of the integral): 

J A,.rdn = j{A, + ^^.r)^\d^x a d^y. 

d\ 

Since A depeds only on x, we can write — A-i, hence the term to be added 

ox'- 



to J A, J' y^ld-^x A d^y is 

J ^^J'y^ld^x AdS^ J {Xr^AG\)■,^ d^x Ad%- J A(J'VN);» A d% 

This term can be written as / div{XJ^)dft — J \div{J^)dft. According to 
the continuity equation and taking into account that A does not depend on y, 
we can write it as J div{XJ^)dft+ J div{XJ^)dft — J div{XJ)dn, i.e., it can 
be written as a boundary term. When performing variations of the action (and 
assuming, as in the classical case, that variations vanish on the boundary), these 
terms will cancel out. 

In conclusion, transformations A{x, y) A{x, y) + dX{x) of the 4-potential 
do not affect the action ([9]). 

Remark. If A; Ai{x), then from ((SH), it follows .P = Q. 



8 Equations of motion 

Let us consider momentarily the case of a single particle. The equations of 
motion are obtained by varying the trajectory x = x{t) in the first two terms 
of ([9]), which are in our case written in the form of a single integral along the 
considered curve: 



S2 = — J {mc^ gij{x, x)x'^x3 -\ — Ak{x, x)x'')dt. (55) 

The 0-homogencity of A insures that the action ^2 is invariant under eventual 
changes of parameter 1 1— > i' of the curve. 
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A further restriction can be imposed on the y-dependence of A in order to 
make all the approach more elegant and provide a simple relation of A with the 
canonical 4-momentum and the equations of motion of charged particles. 

Once the independence of the integral on the parametrization was estab- 
lished, we are free to choose the parameter along the considered curves. Tradi- 
tionally, when deducing the equations of motion, curves are parametrized by the 
arclength s. In this case, the action 5*2 in (|55p is equivalent to the one provided 
by the Lagrangian 

L^^mcg,j{x,y)y'y^ + '^Akix,y)y'', y=^' (^6) 

which is more comfortable in view of Legendre duality and Hamiltonian formal- 
ism. 

The canonical momentum of L is given by 

Pi = = + - {Ak-iy +Ai). 

oy^ c 



In isotropic (pseudo-Riemannian) spaces, if we assume that A ~ A{x), 
then there exists only one potential providing a given interaction Lagrangian 
Lint = Ai{x)y\ But in Finsler spaces, where Ai — Ai{x, y), a given Lagrangian 
Lint = Ai{x, y)y^ can be given by infinitely many functions Ai = Ai{x, y). Thus, 
to a Lagrangian Lint, it corresponds a whole equivalence class of potentials A. 
Comparing to (1191) . it appears as convenient to choose from each class the rep- 
resentative for which 

Afc..y^ - 0. (57) 
We will call this condition upon A, the gradient gauge. In the gradient gauge, 

dy' 

Remark. 0-homogeneity of A insures that we also have Ai.^y^ = 0. 
In the gradient gauge, the canonical 4-momentum is given by 

oy^ c 

dL 

in other words, the Liouville fcanonieal) 1-form 6 — — — rdx"^ attached to L is 

given hy 9 = [mcyi 4- —Ai)dx^ and the Poineare 2-form u = d6, by 
c 

UJ — ^{Aj^i — Anj)dx^ A dx^ — [mcgij -f -Ai.j)dx^ A Sy^ 
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In the following, we assume that the tangent bundle TM is endowed with an 
(arbitrary) Ehrcsmann connection N and the corresponding linear connection 
DT{N), (IMD. 

Variation of (I56p provides the Euler-Lagrange equations: 

ds c ^ c 1 ds' y ds' ^ ' 

Dy^ dy^ 

where — — = —J \- L^^i.y^y'' and we assumed A^.iy^ = 0. 

ds ds ■' 

The first term in the right hand side above is similar to the usual one in 

<? ■ Sy'' 

pseudo-Riemannian spaces, while the second one —F"^^—— is new and appears 

c ^ ds 

due to the dependence of A on the variable y. 

Remark. Both the "traditional" Lorentz force term (given by = -F'^j^y^) 

q 5y-^ 

and the correction given by = -F'^-—r- are orthogonal to the velocity 4- 

c ^ ds 

vector y = X : 

g,,F'y^=Q, g.^F'y^^O. (59) 

Remark. In equations (j58p , we can use any Ehresmann connection N, their 
form does not depend on N. 



9 Stress-energy-momentum tensor 
9.1 In flat pseudo-Finsler spaces 

Let us consider the vector space M — M.* endowed with a flat pseudo-Finsler 
metric 

ffy = ftj(y)- 

Assuming that coordinate transformations are linear (as traditionally done in 
special relativity), we can choose the trivial Ehresmann connection N'^^ ~ 0, 

d d 

hence the adapted frame on TAI is the natural one (Si — — — r, dt — tt-^) and 

ox^ oy^ 

its dual is {dx^,5y^ = dy^). 

Spacetime translations = + i = 0, 3 induce the following transfor- 
mation on TM : 

x'^x' + e\ y'^y\ (60) 

By generalized energy-momentum tensor on TM, we understand the Noether 
current given by the invariance to spacetime translations (accordingly, to trans- 
formations ([50]) ) of the action 

Sf-^- [ -r^F * Fd'^x A d'^y, (61) 

./ IbTTC 
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symmetrized by adding a divergence term. 



As pointed out in Section 2, invariance with respect to transformations (j60l) 
of the above means the absence of exphcit dependence on the base coordinates 
a;* of the Lagrangian. 

In order to find the form of the energy-momentum tensor of the electromag- 
netic field in a flat pseudo-Finslerian space, let us first see how Noethcr theorem 
is read in these spaces. 

For an action 



where A — L^\G\ is a Lagrangian density on TM and — q(^i){x, y) are the 
field variables, the Euler-Lagrange equations are: 

d dk d dk -0 feS) 

The absence of explicit dependence on of A means 

dk _ dk dq(i) dk dk 
dx' dq(i) dx' dq(i)^k oq^iy^ 

dk 

(where we understood also summation over /). Substituting — from (j63p and 

grouping terms, we get 

d f dk 



dx^ 




Thus, the invariance of an action on TM under translations on the (flat) 
base space M leads to the appearance of a quantity consisting of two blocks^ 
namely, by symmetrizing (adding divergence terms to) the quantities 

^fc 1 xfcA ^fc 1 dk 



With these notations, 

d d 



Case 1 (J = 0). In order to "guess" the form of the generalized energy- 
momentum tensor for the electromagnetic field, it is advantageous to assume 
for the beginning that J = and apply the above construction to the Lagrangian 
density: 

A = —F * F = —FbcF^^JWV 

167r IGtt ^ ^' ' 
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oAk^i 47r cf^fe-r 47r 

47r 4 47r 

The obtained Noethcr current can be symmetrized by adding divergence 

terms. It can be easily seen that, for J = 0, the term — (F'^^A^ k + F^^^i-k) 

47r 

can be expressed as a divergence term. 

By adding this term to f\, we get T\ = -^{F^''F^k+F^^F,j,~6\FBcF^^) 
or, equivalently, 

T\ = ^{-F'''F,b + \s\FbcF''^'). (65) 

Similarly, — F'-'^Ai k can be expressed as a divergence term; by adding it to 
47r 



T\ = -^F'-'^Ak r , we obtain 



T'^ - --LiT^'fcF.fc. (66) 
47r 

Thus, in the case J = 0, we get div{T) = 0, i.e.. 

This suggests the following 

Definition 3 The generalized energy-momentum tensor in the flat Finsler space 
{W^,J-{y)) is the symmetric tensor 

T = T,jdx' ®dx' +Tf^dx'- ®dy^ (67) 

with local components given by h65]) and i66\} . 



The horizontal component Tijdx^ ig) dx^ is the the usual energy-momentum 
tensor (plus some correction due to anisotropy), while the mixed one Tfjdx'^®dy^ 
is new. As we have seen above, these new components play a role in the analogue 
of the conservation law. 

Case 2 (J 7^ 0). Let now the TM-current J be arbitrary. By using Maxwell 
equations (with = 0), we get: 

^ '^iT\^\) + -L{T\^\)] = --{F,,J=+F,^J^). (68) 



In brief, 

div{T) ^ ^-ijF. (69) 
c 
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9.2 In general Finsler spaces 

In general (pseudo-)Finsler spaces, we will still define the generalized energy- 
momentum tensor for the electromagnetic field as above: 

T = T.jdx' ®dx^ +T,jdx' ®5y\ (70) 

TrA - ^{-FA''F,B + \9^AFBcF''^), 

An 4 

where g^j — and indices A, B, C take all values corresponding to both hori- 
zontal and vertical components. 

Remark. The horizontal components of the generalized energy-momentum 
tensor can be obtained by varying the action 5*^ with respect to the spacetime 
metric g (i.e., with respect to the horizontal part of the metric (Gab))' 

SgSp ^YcJ T^^t^Sg'^d^l, 

while the mixed components Tq are obtained by varying (independently) Sp 
with respect to the Ehresmann connection N : 



In curved pseudo-Finsler spaces case, the adapted frame {Si,di) is generally 
nonholonomic, hence any linear connection DT{N) which preserves the distribu- 
tions generated by N has generally nonvanishing torsion (at least, Sy^{T{Sk, 5j)) = 
R^jk 7^ 0)- III this case, the covariant divergence of the energy-momentum ten- 
sor is not simply equal to ijF, but has a more complicated expression, 

c 

involving the torsion tensor. The situation formally resembles to the one in 
Riemann-Cartan geometry, [32]. 

In order to find the relation between —ijF and the generalized energy- 

c ^ ^ ' 

momentum tensor, it appears as most comfortable to express — ijF — — (Fij - 

c c 

FfjJ^) in terms of covariant derivatives p4l) . 

Let us assume for simplicity that v = v{x) and N is given by pip . Then, 

Q' = F'''C^^. Q'^'^F^^^R'^k- (71) 
Taking into account the Maxwell equations, we get 

- -^{Fq.P + FqJ^) = + T^.. + T^fi% + T\r\. (72) 
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10 Conclusion 



For a 4-dimensional pseudo-Finsler space (M, F) , we have constructed a notion 
of electromagnetic tensor, based exclusively on variational calculus and exterior 
derivative, [50] . 

The 4-potential is defined as a horizontal f-forni A = Ai{x,y)dx^ on the 
tangent bundle TM, having its components Ai homogeneous of degree in y. 

In terms of this 4-potential, the generalized electromagnetic tensor is the 
2-form F — dA. The Maxwell equations on TM are then written as: 

An 

dF = 0. 5F = J^. 

c 

The TM-current J = FSi + Fdt is a vector field on TM satisfying identically 
div J — 0. Its horizontal component FSt provides the usual notion of 4-current 
(plus a correction term due to the y-dependence of A), while the vertical one 
Fdi is due to the anisotropy of the 4-potential and to the nonholonomy of the 
frame. 

Further, for flat pseudo-Finsler spaces {AI,F{y)), the generalized energy- 
momentum tensor is defined as the symmetrized Noether current corresponding 
to invariance to spacetime translations of the field Lagrangian. We obtained 

T = T.jdx' ® dx^ + Tfjdx' ® dy^ , (73) 

T^A = ^{-F^^Fb + -/aFbcF''^), (74) 

(where (5^ is the Kronecker delta and 5\j — Q and A, B, C take all values corre- 
sponding to both horizontal and vertical components). The generalized energy- 
momentum tensor satisfies the conservation law 

div{T) = --{Fi^J^ + Ff^,P). 
c ■' 

In curved Finsler spaces, the same expressions can be obtained by varying 

by varying the action Sp — — J F * Fd^l for the field with respect to the 

metric gij{x, y) (thus getting Ty) and with respect to the Ehresmann connection 
TV (which provides the components, T!jj). 

The above considerations hold true in more general spaces such as Lagrange 
or generalized Lagrange spaces, with the only mention that in these spaces, the 
equations of motion of charged particles become more complicated (and the 
0-homogeneity assumption on A is dropped). 

Remark. Prior to this model, to our knowledge, there existed only one geo- 
metric model for electromagnetism in Finsler spaces, belonging to R. Miron and 
collaborators, [15], [18], [T7], [19], in which it is defined, by means of deflection 
tensors of linear connections on TM , a notion of electromagnetic tensor on TAI 
(with horizontal hh- and vertical vv- components) and it is provided a general- 
ization of Maxwell equations. Several of the the advantages of our approach are: 



26 



obtaining by variational methods the TM- versions of: Maxweh eqiiations, equa- 
tions of motion (and Lorentz force, respectively) , energy-momentum tensor; an 
easier interpretation of the new (mixed hv-) component of the electromagnetic 
2-form (as appearing, for instance, in the equations of motion); obtaining an 
analogue of the usual continuity equation as an identity, and also, the possibility 
of a compact writing using exterior derivatives. 
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